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Single-point measurements of fluctuations in the scrape-off layer of magnetized plasmas 
are generally found to be dominated by large-amplitude bursts which are associated with 
radial motion of blob-like structures. A stochastic model for these fluctuations is presented, 
with the plasma density given by a random sequence of bursts with a fixed wave form. Under 
very general conditions, this model predicts a parabolic relation between the skewness and 
kurtosis moments of the plasma fluctuations. In the case of exponentially distributed burst 
amplitudes and waiting times, the probability density function for the fluctuation amplitudes 
is shown to be a Gamma distribution with the scale parameter given by the average burst 
amplitude and the shape parameter given by the ratio of the burst duration and waiting 
times. 

Cross-field transport of particles and heat in the scrape-off layer (SOL) of non-uniformly mag- 
netized plasmas is caused by radial motion of blob-like structures [ll|-l6|]. This results in single- 
point recordings dominated by large-amplitude bursts, which have an asymmetric wave form with 
a fast rise and a slow decay, and positively skewed and flattened amplitude probability density 



functions [|7 



14ll . Measurements on a number of tokamak experiments have demonstrated that 



as the empirical discharge density is approached, the radial SOL particle density profile becomes 
broader and plasma-wall interactions increase fllll - llSll . 

Probe measurements on Tokamak a Configuration Variable (TC V) have demonstrated a remark- 
able degree of universality of the plasma fluctuations in the far SOL region, which is dominated 



by radial motion of filament structures and a relative fluctuation level of order unity dlO 



m. In 



particular, the amplitude distribution of the plasma fluctuations are found to be well described by 
a Gamma distribution across a broad range of plasma parameters and for all radial positions in the 
SOL [|12|1 . Excellent agreement was found when comparing analysis of these data with turbulence 
simulations based on interchange motions [l9l-lllll. 

In this Letter, a stochastic model for intermittent fluctuations in the plasma SOL is presented 



with all statistical properties in agreement with experimental measurements. It is demonstrated 
that this model explains many of the salient experimental findings and empirical scaling relations, 
including broad plasma profiles and large fluctuation levels, skewed and flattened amplitude prob- 
ability distribution functions and a parabolic relation between the skewness and kurtosis moments. 
The latter has been observed in the boundary region of numerous experiments on magnetized plas- 
mas as well as in hydrodynamical and astrophysical systems dominated by intermittent fluctuations 



m 



mm. 



There have been several previous attempts at describing the universal features of intermittent 



fluctuations at the boundary of magnetically confined plasmas [|19l- l21|l . However, none of these 
models provide the appealing simplicity, physical insight, novel predictions and favourable com- 
parison with experimental measurements as the theory presented here. In particular, the statisti- 
cal properties which the present model is based on have been directly confirmed by experiments 



m 



m 



Experimental measurements as well as numerical simulations suggest that plasma fluctuations 
in the far SOL can be represented as a random sequence of bursts events, 

^{t)=Y,Akiir{t-tk) (1) 

k 

where A^ is the amplitude and t^ is the arrival time for burst event k, and i/a is a fixed burst wave 
form. This stochastic process resembles a general class of models known as "shot noise", in which 



the noise is generated by the addition of a large number of disturbances [I22N24I1 . The objective 
is to estimate the mean value and higher order moments of 4>, the amplitude probability density 
function P(p and discuss how the burst statistics are related to broad SOL plasma profiles and large 
fluctuation levels. 



If there are K burst events in a time interval T, the average 



T /K. It follows that the mean value of the plasma density is ll22l - l24ll . 



3urst waiting time Tw is given by 



m=^^ r dtw{t). (2) 

Here and in the following, angular brackets are defined as an average of a random variable over all 
its values. The above equation shows that the mean plasma density is given by the average burst 
amplitude and the ratio of the burst duration and waiting times. Equation (O thus elucidates the 
role of burst statistics for high plasma density in the far SOL. It should be noted that this result 
only depends on the integrated burst wave form and the average burst amplitude and waiting time. 



Expressions for the variance and higher order moments of 4> have been derived in the case that 
burst events occur in accordance to a Poisson process with rate I/Tw The probability of exactly 
K burst events in a time interval T is then given by the Poisson distribution, 



P{K) = exp 



w / \ ^w 



K 



1 



(3) 



From this it follows that the 
perimental measurements [|8, 



Durst waiting times are exponentially distributed, as found from ex- 



^t(t) = — exp ( - 



T 



(4) 



The general result states that for the stochastic process defined by Eq. ([T]), the cumulants Kn for the 
probability density /^l^) are given by 



[|23l 



gj 



K 



{A")In 






(5) 



where the integral of the n-th power of the wave form is defined by 



/„=/ dt[xir{t)]'\ 



(6) 



The cumulants are the coefficients in the expansion of the logarithm of the characteristic function 
for P(p, 



ln(exp(z4>M)) = ^ fc^ 



n=l 



{iuY 



(7) 



A power series expansion shows that the characteristic function is related to the raw moments of 
4>, defined by /i^ = (4>"), 



(expo-4>.)) = i+f;i^ = i + f;M:- 



lU] 



n=\ 



n\ 



(8) 



n=l 



Further expanding the logarithmic function in Eq. (|7) and using Eq. ([8]), it follows that the lowest 
order centred moments /i„ = ((4> — (^))") are related to the cumulants by the relations /i2 = K2, 
111, = K2 and jU4 = K4 + 3k-|. 

The variance and higher order moments are straight forward to calculate from Eq. (|5]) for gen- 
eral burst wave forms and amplitude distributions. The coefficient of variation, skewness and 



flatness are given respectively by 
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(9a) 
(9b) 

(9c) 



The two latter relations imply that there is a parabolic relation between the skewness and flatness 
moments, 

(10) 



^^3^/2^(A^,2. 



Such a parabolic relation between the third and fourth order moments have been found for a wide 



m2i\\. 



variety of physical systems dominated by intermittent fluctuations 1121 

The expressions for the higher order moments become particularly simple for a burst wave form 
given by a sharp rise followed by a slow exponential decay, 

(11) 



V/(r) = 0(r)exp 

where is the step function and Td is the burst duration time. This is the typical wave form found 
from probe and gas puff imaging measurements in the far SOL [|7l- ll3|] . The integral given in 
Eq. ((61) is then /„ = Td/n and the cumulants are thus given by K"„ = Td(^")/"'^w- The expressions 
for the coefficient of variation, skewness and flatness become 

.2Td; (A) 

'8Tw\^/2 <A3) 



c 

s 



9rJ (A2)3/2^ 



F = 3+ .. 

The relation between the skewness and flatness is in this case given by 

9(a2)(a4) , 
8 (A3)2 



(12a) 
(12b) 
(12c) 



(13) 



Note that independent of the burst wave form, the probability distribution function for 4> is posi- 
tively skewed, 5 > 0, and flattened, F > 3, for positive definite burst amplitudes A. 



The above expressions for the lowest order moments simplify further in the case of exponen- 
tially distributed burst amplitudes, 



Pa{A) 



1 

(A) 



exp 



{A) 



(14) 



which is also consistent with experimental measurements in the SOL of magnetically confined 

fin 

plasmas [l8i[16|]. The raw amplitude moments are then given by (A") = {A)"n\. In this case the 



relative fluctuation level, skewness and flatness can be written as 
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1/2 



3 + 
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(15) 



'^d / "^d 

All these moments increase with the ratio Tw/Td. The parameter 7= Td/Tw is thus a measure of 
intermittency in the shot noise process. The relation between the skewness and flatness moments 
now becomes 



F = 3 + -S^ 



(16) 



which is in excellent agreement with measurements in the SOL of tokamak plasmas |12L II6L 



m. 



The interpretation of these results is evident. For short waiting times and long burst duration, 
the signal 4> will at any time be influenced by many individual bursts, resulting in a large mean 
value and small relative variation. In the opposite limit of long waiting times and short duration, 
the signal is dominated by isolated burst events, resulting in a smaller mean value and large relative 
fluctuations, skewness and flatness. This is clearly illustrated in Fig.[Tl which shows two numerical 
examples of the shot noise process given by Eq. ^ for exponentially distributed burst waiting 
times and amplitudes. 

The results presented above show that the skewness and flatness vanish in the limit of large 
7. It can be demonstrated that the probability density function for $ then approaches a normal 



distribution [123 , 
given in Eq. d?]). 



2411 . The distribution Pep can be written in terms of the characteristic function 



P'^m 
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27r 



du exp 



?4>a + ^ 



n=l 



(17) 



In the limit of small T^/Td the exponential function can be expanded as a power series in u. 
Integrating term by term then gives 



^rms/'*(4>) 



1 



(2;r)i/2 



exp 




1 + 



M3 



3!4>3„,(2;r)i/2 



(4)3-34))+i?($) 



(18) 



where R is the sum of the remaining terms in the expansion and the centred and rescaled amplitude 
is defined by ^ = (4> — (4>))/4>]-nis- The terms inside the square bracket in Eq. ([18} are of order 
1, (Tw/'T'd)^'^^ and Tw/Td, respectively. This result shows how the probability density function for 
4> approaches the normal distribution in the limit of large 7. This transition to normal distributed 
fluctuations is expected from the central limit theorem, since in this case a large number of burst 
events contributes to 4> at any given time. 

The asymptotic probability density function in the strong intermittency regime for small 7 can 
be obtained by neglecting overlap of individual burst events. Considering first a single burst event 
(j){t) = Aexp(— r/id), the time dt spent between (j) and + J(^ is given by dt/d^ = Td/(j). Note that 
due to the assumed exponential wave form, the burst amplitude A does not enter this expression. 
The number of bursts with amplitude above 4> is given by the complimentary cumulative ampli- 
tude distribution function, which for an exponential distribution is exp(— $/(A)). The probability 
density function P^ is given by the proportion of time which $(f ) spends in the range from 4> to 
4> + (i4>. With the appropriate normalization, the asymptotic probability density function in the 
strong intermittency regime is thus given by 

limPc^(4>) = lim-i-;^exp(--^), (19) 

where we have defined the Gamma-function 

/•CO 

r(7)= / J(p(p^-iexp(-(p), (20) 

Jo 

which to lowest order is given by 1 /7in the limit of small 7. This probability density function has 
an exponential tail for large amplitudes but is inversely proportional to $ for small amplitudes due 
to the long quite period between burst events in this strong intermittency regime. 

The characteristic function for a sum of independent random variables is the product of their 
individual characteristic functions. Thus, the probability that a sum of A' burst events ^^ lies in the 
range between 4> and 4> + J4> is given by 



27r 



/CO K. 

jMexp(— ?4>m) J^ (exp(z0^M)), (21) 



k=\ 

where the characteristic functions are averaged over the values of (^k- For general amplitude dis- 
tribution and burst wave forms, 

{oxp(i^ku)) = - dtk dAPAiA) exp [iAuXf/it - tt)] , 
I Jo J-00 



where T is the duration of the time interval under consideration. Since all the K characteristic 
functions in Eq. (|2TI) are the same, the conditional probability P]^ is given by 



1 /.oo 

2k 7-00 



assuming the number of events ^ in a time interval T to be given. The probability density function 
for the amplitude 4> is given by summing over all K, 

oo 

P^{^)= Y,P{K)PKm, (22) 

K=Q 

where P{K) is given by Eq. ([3]) The stationary probability density function for $ is obtained by 
letting T —^ oo. Some elementary manipulations lead to the desired result, 

p^(4)) ^— f exp i -?4>M H / dAPA{A) f dt [exp(zAMV/(f )) - 1] i ■ (23) 

^/t J — oo [ Tw J — °° J — °° J 

The logarithm of the characteristic function of Pt^ is thus 

1 f^ f°° oo 1 (ill)" /■°° f'°° 

— dAPA{A) dt[Qxp{iAu\ir{t))-l\=Y,—^4- dAA"PA{A) dt[\ir{t)]". (24) 

where again the exponential function has been expanded. This establishes the general result stated 
by Eq. ©. 

For the special case of exponentially distributed burst amplitudes, the amplitude integral in the 
above equation is given by {A)"nl, cancelling the factorial in Eq. (I24l ). Further invoking the ex- 
ponential wave form given in Eq. (ITT)) , it follows that the characteristic function for the stationary 
distribution can be written as 



exp 



^,f {i{A)u) 



n=\ 



{\-i{A)uy^. (25) 



This is nothing but the characteristic function for a Gamma distribution with scale parameter (A) 
and shape parameter 7. Thus, the probability density function for 4> is given by 

The lowest order moments and asymptotic limits of this distribution agree with the expressions 
discussed previously. In particular, the mean value is given by (4>) = 7(A) and the variance by 



8 

For 7 > 1, the most likely amplitude of the Gamma distribution is (7— 1)(A) and the shape 
of the distribution function is unimodal and skewed. When 7= 1, Pep becomes an exponential 
distribution with the mean density given by the average burst amplitude, 

''*(*) =(i)«p(- (I))- <''' 

Note that by writing the average burst amplitude as (A) = (^)/7, the Gamma distribution given in 
Eq. (|26|) can be written in terms of the average plasma density as 

($)P^(4.) = -I- f ^y expf-^y (28) 



r(7)VW; "V w 

where the scale parameter is given by (4>)/7and the shape parameter is 7= (4>) /^^s- ^^ should 
be noted that there are no fit parameter when comparing this prediction to experimental measure- 
ments. The above equation is exactly the form of the Gamma distribution found empirically to 
describe plasma fluctuations in the SOL of TCV across a broad range of plasma parameters [I12I1 . 
Recently, it has also been shown to describe fluctuations recorded by gas puff imaging measure- 
ments in the SOL of the Alcator C-Mod tokamak [lldl . 

In summary, a stochastic model for intermittent fluctuations in the boundary region of magne- 
tized plasmas has been constructed as a random sequence of bursts which represent radial motion 
of blob-like structures. The mean plasma density is given by the average burst amplitude and the 
ratio of burst duration and waiting times. In the case of exponentially distributed burst amplitudes 
and waiting times, the amplitude probability density function is shown to by a Gamma distribution. 
This simple model thus explains the salient fluctuation statistics found in numerous experimental 
measurements and elucidates the role of burst statistics for large SOL plasma densities and fluctu- 
ation levels. The general parabolic relation between skewness and kurtosis moments predicted by 
this model likely explains the wide spread observation of this scaling relation in physical systems 
dominated by intermittent fluctuations. 

Discussions with B. LaBombard, M. Melzani, H. L. Pecseli, R. A. Pitts, M. Rypdal and 
J. L. Terry are gratefully acknowledged. 
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FIG. 1. Simulated shot noise time series for exponentially distributed burst amplitudes and waiting times. 



